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Abstract 

We calculate the ground state properties of recently synthesized superheavy nuclei starting from Z=105- 

120. The nonrelativistic and relativistic mean field formalisms is used to evaluate the binding energy, 

charge radius, quadrupole deformation parameter and the density distribution of nucleons. We analyzed the 

stability of the nuclei based on the binding energy and neutron to proton ratio. We also studied the bubble 

^ structure of the nucleus which reveals about the special features of the superheavy nucleus. 
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PACS numbers: 



21.10.Gv, 21.10.Tg, 21.10.Dr, 21.60.-n 



I. INTRODUCTION 

The study of superheavy nuclei is an interesting topic for the current days research in nuclear 
physics. Initially, the transuranium elements were made by subjecting Uranium to a high neutron 
flux allowing to capture neutrons successively. The neutron-rich isotopes undergo /3— decay and 
produce new elements. This process allows to prepare element maximum upto Fermion {Z = 
100), because of the shorter time interval for the spontaneous fission of the new element than 
neutron capture time. To make heavier Z nuclei one has to reach the so called island ofstability 
which is a matter of discussion from last five decades II -4]. Earlier, it was predicted ll], [sl that 



the next proton magic number would be Z=\26, because of the equality of the proton and neutron 



magic numbers for known closed shell nuclei. However, several microscopic calculations 



suggest a shift of this number toll4. One of the cause of the shift is Coulomb effects on the 



spherical single particle levels. The use of shell correction by V. M. Strutinsky 11211 to the liquid- 
drop calculation of binding energy (BE) opens a more satisfactory exploration of double close 
nuclei beyond ^°^Pb. Recently, based on the relativistic mean field (RMF) and Skyrme Hartree- 



Fock (SHE) calculations Patra and collaborator IIS 



IM predicts that the island ofstability in the 
superheavy region is situated near mass number A ~ 300 peaking at ^^^'^^"^120. 

Now the main issue is how to reach the island of stability ? Because of the limitation of the 
nucleosynthesis by neutron capture, the only possible way to reach the stability island is nuclear 
reaction induced by heavy ions via complete fusion reactions the so-called coldfusion reactions 



lilSh . In this method, the kinetic energy of the heavy projectile is absorbed in a compound nu- 



cleus (heavy target). It is interesting to note that the superheavy elements with Z=105-118 are 



synthesized using cold-fusion reactions lild llTII. Although, the superheavy elements Z=93-118 



are synthesized artificially, the knowledge about their chemistry is limited to Z < 107 due to the 



shorter /3— decay life time Tjs UM- The chemistry of elements Z=101 and higher Z nuclei were 
unknown and the prospects of the studies of their properties did not look promising because of 
their shorter life time. However, some isotopes of these elements are recently discovered whose 
life time seems to be relatively large, which gives a hope for the studies of the chemistry of the su- 
perheavy elements. In this paper, our aim is to study the properties and the stability of the recently 
synthesized superheavy nuclei. 



The stability of a nucleus mostly depends on its structure which is again a function of the 
combination of number of neutron and proton. For an extremely large nucleus as it is in super- 
heavy region, the shape of the ground state configuration most unlikely spherical. Because of the 
large Coulombic repulsion, the central portion of the nucleus may take a depletion and create a 
bubble like structure. There are some interesting and extensive calculations for various proper- 
ties including bubbles or semi-bubbles structure of superheavy elements [1191 120I1. In this paper, 
the stability of the superheavy nuclei are discussed by two well known self-consistent mean field 
models, namely Skyrme Hartree-Fock (SHF) using SkI4 parameters [12 ill with BCS pairing and 
Relativistic Mean Field (RMF) model with NL3* parameter set 12211 . 

The paper starts with a short introduction in Sec. I. The formalisms of Skyrme Hartree-Fock 
and relativistic mean field theory are presented in Sec. II. Mostly, the stability and structure of the 
synthesized nuclei will be discussed in the Results and discussions are given in Sec. III. In this 
section, mostly the stability and structure of the synthesized nuclei will be discussed. The bubble 
properties will also be analyzed taking into consideration the density distributions of proton and 
neutron. Sec. IV contains a brief summary of the paper. 



n. THEORETICAL FORMALISMS 



A. Skyrme Hartree-Fock (SHF) method 



There are many known parametrization of Skyrme interaction which reproduce the experimen- 
tal data for ground state properties of finite nuclei |123i |240 as well as the properties of infinite 
nuclear matter upto high density \2^. The general form of the Skyrme effective interaction can be 



expressed as a density functional "H with some empirical parameters ^ 



2],y, 



H — K, -\- + TH-z + 'Heff + 



(1) 



where /C is the kinetic energy. Ho the zero range, Hz the density dependent and Heff the effective- 
mass dependent terms, which are relevant for calculating the properties of nuclear matter. More 
details can be found in Refs. [|24l.l26l. l27|l . These are functions of 9 parameters U, Xi {i = 0, 1, 2, 3) 



and 77 are given as 



Uo = -^to [(2 + xo)p' - (2x0 + l)(Pp + pD] , 
^3 = ^hp' [(2 + xs)p' - (2x3 + l)(pj + pI 



Heff = -[tl(2 + Xi)+t2(2 + X2)]rp 
o 

+ 1 [^2(2X2 + 1) - tl(2Xi + 1)] {TpPp + TnPn). 



fc2 

The kinetic energy /C = y-r, a form used in the Fermi gas model for non-interacting Fermions. 



(2) 
(3) 



(4) 



2m 



The other terms T-Lsp and H^j with 64 and 64 are the surface contributions which also responsible 
for a better spin-orbit interaction are defined as 



n 



Sp 



1 

16 



3h{l + -x,] 



h{l + -X2) 



(Vp)^ 



1 

16 



1 , 1. 
3ti(xi + 2) + ''^2(3^2 + -, 



{VPnf + (Vp, 



and 



SJ 



hipV ■ J + h'^{pnV ■ Jn + PpV ■ Jp 



(5) 
(6) 



Here, the total nucleon number density p = Pn + Pp, the kinetic energy density r = r„ + rp, and the 
spin-orbit density J = Jn + Jp. The subscripts n and p refer to neutron and proton, respectively, 
and M is the nucleon mass. The Jg = 0, q = n or p, for spin- saturated nuclei, i.e., for nuclei with 
major oscillator shells completely filled. The total binding energy (BE) of a nucleus is the integral 
of the density functional H. In our calculations, we have used the Skyrme Ski4 parameter set [12 ill . 
This set is known to well account the spin-orbit interaction in finite nuclei, related to the isotope 
shifts in Pb region and is better suited for the study of exotic nuclei. Several more recent Skyrme 
parameters such as SLyl-10, SkX, SkI5 and SkI6 are obtained by fitting the Hartree-Fock (HF) 
results with experimental data for nuclei starting from the valley of stability to neutron and proton 



drip-lines 



27 



B. Relativistic mean field (RMF) formalism 



From last few decades, the RMF theory is applied successfully to study the structural properties 
of nuclei throughout the periodic table ll30l - l33ll starting from proton to neutron drip-lines. The 



starting point of the RMF theory is the basic Lagrangian containing nucleons interacting with a—, 
cu— and p— meson fields. The photon field is included to take care of the Coulomb interaction 
of protons. The relativistic mean field Lagrangian density is expressed as 130-34]. 



- Qpi^a^fiJiR^^ - ei>iY ^^~J^'\ iA^. (7) 

Here M, m^, rrii^ and nip are the masses for nucleon, a-, to- and p-mesons and ip is its Dirac 
spinor. The field for the a-meson is denoted by a, w-meson by and p-meson by R^. Qs, g^, Qp 
and e^/47r=l/137 are the coupling constants for the a, u, p-mesons and photon respectively. g2 
and (73 are the self-interaction coupling constants for a mesons. By using the classical variational 
principle we obtain the field equations for the nucleons and mesons. 

{- A +m^}(7°(r^,z) = -gapsir±,z) 

- 92(^'^ir^,z) - g3a^{r±,z), (8) 

{- A+m^}l^°(r^,z) = gu.p.{r^,z), (9) 

{- A +mj}i?°(r^, z) = gppsir^, z), (10) 

~AA\r^,z) = ep,{r^,z). (11) 

The Dirac equation for the nucleons is written by 

{-ta.y + V{r^,z)+m^}tP^ = ei^i. (12) 
The effective mass of the nucleon is 

= M + S{r^, z)=M + gMru z), (13) 

and the vector potential is 

Vir^, z) = g^V\r^, z) + gp^R^r^, z) + e^-^p^A\r^, z). (14) 

A static solution is obtained from the equations of motion to describe the ground state properties of 
nuclei. The set of nonlinear coupled equations are solved self-consistently in an axially deformed 



hannonic oscillator basis Np = Nb = 20. The quadrupole deformation parameter is extracted 
from the calculated quadrupole moments of neutrons and protons through 



Q Qn ~l~ Qf 

where R= l.2A^/^. 

The total energy of the system is given by 



167r / 3 
47r 



AR^f32 



(15) 



total 



Epart + Ea- + E^ + Ep + Ec + Epair + E, 



(16) 



where Epart is the sum of the single particle energies of the nucleons and E^, E^, Ep, Ec, Epair, 
Ecm are the contributions of the meson fields, the Coulomb field, pairing energy and the center- 
of-mass energy, respectively. We use the recently reported NL3* parameter set [|22n in our calcu- 
lations for RMF formalism [221. 



III. RESULTS AND DISCUSSIONS 



In the present paper, we calculate the binding energy, root mean square (rms) charge radius 
rch and quadrupole deformation parameter (32. A special attention is given to the neutron and 
proton density distributions of these superheavy nuclei, where we search for the bubble structure 
at various excited intrinsic states. The results are explained in subsection III.A to III.E. 



A. Stability of nuclei 

While analyzing the nuclear landscape, the stability of nucleus is maximum with a neutron 
to proton ratio N/Z at about 1 for light mass region. This value of N/Z increases with mass 
number A to neutralize the Coulomb repulsion and saturated at ~ 1.5. For example, the neutron to 
proton ratio for ^"'^Pb is 1.537 and for 235,238y values are 1.554 and 1.587 respectively. This 
saturation nature of N/Z ~ 1.5 suggests the formation of superheavy nuclei with an N : Z =1.5. 
To see the stability of the synthesized isotopes in the superheavy region, we plot the N/Z values 
in Figure 1 for Z=\05 to 1 18 in addition to the predicted double close nuclei 292,293,304]^ 20. To our 
surprise, we find the minimum and maximum N : Z as 1.433 and 1.614, respectively. 

Another important quantity is the binding energy (BE) of a nucleus. This is the only experi- 
mental observable which can be measured precisely and tells many things about the stability. A 
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FIG. 1: The N/Z ratio of the recently synthesized nuclei with atomic number Z = 105 — 118. The N/Z 
values of the predicted superheavy isotopes 292,293,304^20 are also given. 

nucleus is found to be more stable with a larger binding energy per particle {BE/ A). The trend of 
BE/ A as a function of mass number A is shown in Figure 2. It is clear that low mass nuclei in the 
valley of /3— stability have more life-time compared to the heavier isotopes. This may be due to the 
large Coulomb repulsion of the superheavy nuclei, although the N/Z ratio is almost constant. The 
declining trend of BE/ A indicates the difficulty in the formation of heavier superheavy nuclei. Of 
course, it is to be recalled that the formation of a neutron-proton bound state system requires only 
2.224 MeV which is the typical binding energy of Deuteron {BE /A ~ 1.112 MeV) Bsll . 



B. Shape and shape co-existence 



Nuclei have different binding energy corresponding to different shape configuration in the 
ground or excited states. In some cases, it so happens that the binding energy of two different 
configurations coincide or very close to each other known as shape co-existence |38-4l|. This 
phenomenon is more common in superheavy region, leads to complex structure of these nuclei 



TABLE I: The BE (MeV), /32, Vch (fm) (ground state) and D.R (%) using SHF(SkI4). 



Nucl. 
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BE 
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/52 


rch. 


Spherical 


Prolate 


Oblate 


FRDM 




Sphe. 


Prolate 


Oblate 










D.R„ 


D.F.p 


D.R„ 


D.Rp 


D.R„ 


D.Fp 


BE 


/32 


266Db 




1950.3 


1939.7 




0.262 


-0.282 


6.08 


0.0 


0.0 


0.0 


2.98 


1.32 


7.86 


1953.42 


0.221 






1956.5 


1945.9 




0.258 


-0.283 


6.09 


0.0 


3.41 


0.0 


3.41 


1.29 


7.43 


1960.49 


0.230 


268Db 




1961.7 


1951.3 




0.254 


-0.193 


6.09 


0.0 


3.42 


0.0 


3.83 


1.35 


7.36 


1965.86 


0.221 


269Db 




1966.9 


1957.6 




0.250 


-0.191 


6.10 


0.0 


3.58 


0.0 


3.60 


1.13 


7.53 


1972.01 


0.221 


270Db 




1971.9 


1963.7 




0.259 


-0.188 


6.11 


0.0 


3.65 


0.0 


2.12 


0.81 


7.64 


1976.81 


0.212 


258sg 




1890.5 


1880.0 




0.287 


-0.287 


6.05 


2.73 


4.52 


2.73 


4.54 


2.74 


0.89 


1895.22 


0.248 


260sg 




1905.5 


1894.9 




0.284 


-0.286 


6.06 


2.29 


4.24 


0.0 


4.24 


1.75 


12.86 


1909.90 


0.239 


262sg 




1920.1 


1909.3 




0.279 


-0.283 


6.07 


1.05 


3.69 


1.72 


3.70 


1.59 


12.19 


1923.94 


0.229 


^^^Sg 




1976.5 


1968.0 




0.254 


-0.188 


6.12 


1.39 


4.12 


0.0 


3.45 


0.93 


7.71 


1981.87 


0.212 






1969.8 


1958.7 




0.254 


-0.192 


6.11 


0.80 


4.00 


0.0 


4.02 


1.53 


7.57 


1973.53 


0.222 






1980.8 


1971.9 




0.244 


-0.187 


6.12 


0.0 


4.21 


0.0 


4.62 


1.07 


7.82 


1985.27 


0.221 






1987.4 


1984.7 




0.237 


-0.182 


6.13 


0.0 


4.25 


0.0 


3.16 


0.52 


8.20 


1996.40 


0.192 


264Hs 




1923.6 


1912.8 




0.271 


-0.283 


6.08 


0.0 


4.42 


1.23 


3.74 


1.36 


12.27 


1927.62 


0.229 


268Hs 
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1941.3 




0.262 


-0.281 


6.12 


0.0 


3.60 


0.0 


4.40 


1.26 


11.46 


1957.28 


0.230 


270Hs 




1967.4 


1955.2 




0.256 


-0.193 


6.11 


0.0 


4.17 


0.0 


4.17 


1.76 


7.69 


1971.42 


0.231 






1978.9 


1968.8 




0.248 


-0.188 


6.12 


0.0 


4.40 


0.0 


4.40 


1.48 


7.84 


1984.20 


0.222 


275Hs 




1995.8 


1988.7 




0.233 


-0.180 


6.14 


0.0 


4.39 


0.0 


3.84 


0.78 


8.42 


2000.78 


0.183 


274Mt 




1987.4 


1978.7 




0.237 


-0.183 


6.14 


0.0 


4.71 


0.0 


4.72 


1.49 


8.28 


1991.77 


0.222 


275Mt 




1993.3 


1985.5 




0.232 


-0.181 


6.14 


0.0 


4.31 


0.0 


4.73 


1.27 


8.49 


1998.36 


0.212 






1999.1 


1992.2 




0.225 


-0.178 


6.14 


0.0 


4.27 


0.0 


4.67 


1.05 


8.73 


2003.64 


0.202 


278Mt 




2010.6 


2005.5 




0.212 


-0.173 


6.15 


0.0 


4.38 


0.0 


4.49 


0.70 


9.30 


2015.32 


0.136 


270Ds 




1955.7 


1945.5 




0.250 


-0.195 


6.12 


0.0 


4.30 


0.0 


4.31 


2.09 


8.13 


1960.53 


0.221 


279Ds 


2014.0 


2011.7 


2009.1 


0.199 


0.276 


-0.171 


6.20 


0.0 


4.29 


0.0 


4.81 


1.08 


9.70 


2019.39 


0.127 


28lDs 


2019.8 


2021.8 


2022.3 


-0.030 


0.340 


-0.166 


6.24 


16.96 


22.10 


0.0 


4.44 


0.88 


10.34 


2031.91 


0.108 


277Rg 




1997.7 


1991.3 




0.220 


-0.176 


6.16 


0.0 


4.78 


0.0 


4.32 


2.03 


9.38 


2003.04 


0.193 


278Rg 




2004.0 


1998.4 




0.220 


-0.173 


6.17 


0.0 


4.12 


0.0 


4.72 


1.86 


9.63 


2008.68 


0.164 


279Rg 




2007.1 


2005.4 




0.288 


-0.170 


6.21 


0.0 


4.09 


0.8 


4.24 


1.68 


9.90 


2015.73 


0.164 
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Continue of table-I 
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BE 
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/32 


^2 


Teh. 


Spherical 


Prolate 


Oblate 


FRDM 




Sphe. 


Prolate 


Oblate 










D.F.„ 


D.F.p 


D.F., 


D.F.p 


D.F„ 


D.Fp 


BE 


/32 


280Rg 




2012.5 


2012.3 




0.334 


-0.167 


6.24 


0.0 


4.11 


0.0 


4.83 


1.53 


10.19 


2021.81 


0.117 


28lRg 


2016.9 


2018.4 


2019.2 


-0.036 


0.353 


-0.165 


6.26 


17.27 


22.44 


0.77 


4.56 


1.40 


10.50 


2028.85 


0.099 


282Rg 


2023.8 


2024.3 


2025.9 


-0.037 


0.366 


-0.163 


6.27 


17.06 


20.14 


0.62 


4.61 


1.29 


10.79 


2034.82 


0.099 


282112 


2020.4 


2020.4 


2022.3 


-0.030 


0.370 


-0.161 


6.28 


17.85 


23.20 


0.0 


4.75 


1.94 


11.07 


2032.68 


0.089 


283112 


2027.9 


2028.2 


2029.2 


-0.040 


0.489 


-0.159 


6.38 


16.73 


22.14 


0.0 


3.99 


1.29 


11.33 


2038.67 


0.089 


28M12 


2035.1 


2035.2 


2036.0 


-0.066 


0.525 


-0.159 


6.41 


14.18 


20.42 


1.41 


2.78 


2.34 


11.93 


2045.88 


0.089 


285112 


2042.1 


2041.6 


2042.5 


-0.068 


0.532 


-0.167 


6.43 


13.29 


19.93 


0.0 


2.47 


1.16 


11.20 


2051.61 


0.089 


29^112 


2099.4 


2093.3 


2098.9 


-0.004 


0.559 


-0.005 


6.50 


3.07 


15.17 


0.0 


5.74 


0.22 


11.86 


2103.92 


0.00 


282113 


2021.0 


2021.1 


2017.9 


-0.164 


0.171 


-0.159 


6.18 


1.02 


9.57 


0.0 


8.43 


2.71 


11.44 


2027.20 


0.071 


284113 


2031.0 


2034.4 


2032.3 


-0.002 


0.158 


-0.154 


6.18 


18.47 


24.09 


0.0 


8.61 


2.37 


11.93 


2040.95 


0.080 


285113 


2038.2 


2041.1 


2039.2 


-0.023 


0.151 


-0.155 


6.18 


18.02 


23.72 


0.0 


8.79 


2.06 


11.90 


2048.18 


0.072 


286113 


2045.3 


2047.7 


2045.9 


-0.021 


0.144 


-0.163 


6.19 


17.54 


23.44 


0.0 


8.94 


1.58 


11.65 


2054.22 


0.072 


286114 


2041.4 


2041.1 


2042.0 


-0.012 


0.555 


-0.150 


6.46 


18.59 


24.53 


0.0 


0.0 


2.63 


12.55 


2051.59 


-0.096 


287114 


2048.6 


2047.7 


2048.8 


-0.012 


0.554 


-0.158 


6.46 


18.06 


24.19 


0.0 


4.82 


2.06 


12.13 


2057.65 


-0.078 


288114 


2055.8 


2054.1 


2055.5 


-0.024 


0.547 


-0.165 


6.46 


16.76 


23.18 


0.0 


4.63 


2.14 


12.09 


2065.01 


0.053 


289114 


2062.8 


2060.3 


2062.2 


-0.011 


0.537 


-0.149 


6.46 


16.07 


22.95 


0.0 


4.81 


2.23 


12.36 


2071.04 


-0.052 


296114 


2108.9 


2102.3 


2108.4 


-0.003 


0.565 


-0.043 


6.53 


4.05 


16.58 


0.0 


4.32 


1.68 


13.68 


2113.22 


0.00 


298114 


2121.1 


2113.3 


2120.8 


0.000 


0.566 


-0.012 


6.54 


0.0 


11.18 


0.0 


4.72 


0.0 


10.74 


2123.30 


0.00 


287115 


2043.1 


2044.9 


2044.1 


-0.008 


0.153 


-0.137 


6.20 


19.86 


25.48 


0.0 


6.57 


4.18 


13.69 


2052.72 


-0.096 


288115 


2050.5 


2051.9 


2051.2 


-0.005 


0.145 


-0.137 


6.20 


19.24 


25.09 


0.0 


7.07 


3.66 


13.38 


2059.12 


-0.087 


289115 


2057.8 


2058.8 


2057.9 


-0.002 


0.137 


-0.141 


6.21 


18.19 


25.44 


0.0 


7.48 


3.02 


12.92 


2066.45 


0.053 


290115 


2064.9 


2064.7 


2065.2 


0.000 


0.117 


-0.119 


6.20 


16.88 


23.61 


2.34 


10.75 


3.83 


13.59 


2072.59 


-0.079 


291115 


2071.9 


2071.8 


2072.4 


0.000 


0.064 


-0.102 


6.20 


15.46 


22.70 


8.62 


19.14 


4.39 


14.44 


2079.83 


-0.061 


290116 


2059.5 


2060.5 


2060.3 


0.002 


0.134 


-0.121 


6.21 


17.92 


24.59 


0.0 


7.90 


4.92 


15.51 


2068.76 


0.072 


29^116 


2066.7 


2066.7 


2067.5 


0.002 


0.110 


-0.108 


6.21 


16.61 


23.78 


0.3 


11.81 


5.56 


15.29 


2074.84 


0.072 


292116 


2073.9 


2074.1 


2074.7 


0.002 


0.062 


-0.097 


6.21 


15.32 


22.96 


10.18 


19.89 


5.37 


15.45 


2082.48 


-0.070 


293116 


2080.9 


2081.1 


2081.8 


0.002 


0.052 


-0.092 


6.21 
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FIG. 2: The binding energy per nucleon {BE/ A) is compared with the experimental data. The stars and 
circles are for RMF(NL3*) and SHF(SkI4), respectively. 

which provides the information about the oscillation of nuclei between two or three existing shape 
by perturbing small energy. Here, in our studied nuclei, we find many such examples, where the 
ground and first excited binding energies are degenerated. Also, in several cases the nuclei have 
well defined three distinct solutions (oblate, spherical, prolate) with almost same energy which 
strongly indicates the nuclei may oscillate from oblate to spherical to prolate configuration. 

For analyzing the result , we took a small binding energy difference around < 2 MeV to take 
care of shape co-existence for our calculations. Because of this small difference in binding energies 
between these two solutions, the ground state can be changed to the low-lying excited state or vice- 
versa by a small change in input parameters like pairing gaps etc. It is seen from Tables I and II the 

shape co-existence is observed for 277,278,281,282Rg^ 282,283,284,285 ^12, 286,287,288114 294,297118 ^^d 

292,293120 in SHF calculations and RMF also produces the shape co-existence for ^^^112, ^^^113, 

288,289,296,298114^ 290,291ii5^ 291,29211^^ 293,294,297 1 17^ 297ii8^ 292,293i20. A shapC CO-CXistCnce of 

nuclei indicates that there is a competition among the different shapes to acquire the possible 
ground state for maximum stability and the final shape could be an admixture of these low-lying 
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C. Binding energy and quadrupole deformation parameter 

The ground state binding energy and quadrupole deformation parameter for all the observed 
nuclei considered here are given in Tables I and 11. The much discussed superheavy isotopes 
292,304 j^20 which are considered to be the next double closed nuclei beyond 208pb are also listed. 
To find the ground state solution, the calculations are done with an initial spherical, prolate and 
oblate quadrupole deformation parameter (3q both in the Skyrnie Hartree-Fock (SHF) and relativis- 
tic mean field (RMF) formalisms. The well known SkI4||2ll] and NL3*[|221 parametrizations are 
used through out the SHF and RMF calculations, respectively. It is to be noted that the maximum 
binding energy corresponds to the ground state and all other solutions are the intrinsic excited state 
configurations. Considering this criteria, we evaluate the ground state binding energy and also the 
first and second excited states wherever available. The results obtained from SHF are listed in 
Table I and those for RMF are listed in Table II. Since the experimental binding energy of these 
nuclei are not known, to get a first hand information about the predictive power of our calculations, 
we have inserted the most commonly used binding energy of finite range droplet model (FRDM) 



[2^]- The quadrupole deformation parameter (32 obtained in FRDM is also given for comparison 



3711 . Comparing Tables I and II, it is remarkable to note that the binding energy obtained by 



TABLE II: Same as Table I with RMF(NL3*) model. 
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FIG. 3: The neutron, proton and total matter density distribution for some of the selected superheavy nuclei 
using SHF(Skl4). 

SHF(SkI4) and RMF(NL3*) are in excellent agreement agreement with FRDM prediction. The 
deviation of our results (from the FRDM prediction) is a maximum of about 0.5%. A further in- 
spection of the results finds an edge of RMF(NL3*) over SHF(SkI4). In general, the calculated 
binding energies and deformation parameters are almost agree with each other in all the three 
models. 

D, Density and bubble structure 

The density of a nucleus has the gross information about the size, shape and distribution of nu- 
cleons. The density distribution from SHF(SkI4) and RMF(NL3*) for some of the selected nuclei 
are given in Figures 3 — 7 obtained for all the three solutions, i.e. spherical, prolate and oblate. 
Normally, the density at the center of nucleus has the maximum value and decays to zero at the 
surface. However, in some specific cases, e.g. for spherical solutions of ^^^113, ^^^117, 292,293^^20, 
this trend of density distribution shows an anomalous behavior, i.e. a dip at the center and a hump 
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nearby to it following a slow decreasing in density to zero at the surface. This type of density 
distribution is known as bubble structure [4^. The first possibility for existence of toroidal or 
bubble structure other than the spherical topological structure was sug gest ed by Wheeler [|42ll . 
The occurrence of bubble nuclei has been extensively studied by Wilson [|43|l and later by Siemens 
and Bethe [|44l] . Several models like independent particle model [|45|1 and Hartree Fock Model [|46(1 
also investigated the possibility of low density region at the center of the nuclei. This structure is 



481] 



not confined to a particular region, but have the possibility for light to superheavy nuclei [|47l . 
The contribution of density at r = is offered by non zero wavefunction of 5- states and depop- 
ulation of this level leads to the depletion of central density and formed a bubble like structure 
148(1 . The inversion of 5-states to other higher states may be the possibility of depletion of central 
density and can be interpreted as s-d orbital inversion l|43,0]- 

The bubble effect can be realized quantitatively by calculating depletion fraction for the neu- 
trons as well as protons using the relation ll48n 



= ^^'"""^^ 1^^^"^" X 100, 

\Pmax ) a 



(17) 



where pmax, pcen represent the maximum and central density respectively and a denotes neutron 
or proton. The calculated values of depletion fraction for neutron as well as proton by using non- 
relativistic and relativistic densities are given in Tables I and II for all three solutions. In our 
present study, most of the nuclei achieve the prolate shape as their ground state solution both in 
the SHF and RMF calculations. Here, we do not find significant bubble structure, i.e. no good 
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FIG. 4: Same as Fig. 3, but for 292,293,304^20 with SHF(SkI4). 



amount of depletion fraction are estimated for this shape, except a few cases. It is evident from 
Tables I and II, a remarkable depletion fraction is there in 292,293 ]^20 over the whole series of nuclei 
on spherical solution in both models. 



E. Qa and for superheavy nuclei 

The Qa value of a nucleus provides valuable information for its stability. This is more im- 
portant for superheavy nuclei, which is the central theme of study in the present paper. Thus, it is 
worthwhile to evaluate the Qa values for the nuclei which will also give us to estimate the life time 
Ta of the considered nuclei. It is to be noted that a— decay is the decay mode for these elements 
(and not the /3— decay) which lie on the valley of stability with a neutron to proton ratio ~ 1.5 (See 
Fig. 1). The calculated Qa and Ta for Z = 105 - 120 within the relativistic and non-relativistic 
formalism using NL3* and SkI4 parameters are listed in Table III. The Qa is calculated by using 
the following expression [Isill 

Qa = BE{N, Z) - BE{N -2,Z -2)- BE{2, 2). (18) 



TABLE ni: The calculated and Ufe time using SHF(SkI4) and RMF(NL3*) are compared with 
FRDM and experimental data. 
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9.04 


9.61 




14.42s 


2.63s 


10-^-26, 






8.84 


8.88 


8.40 


8 54±0 08 


125 06s 


95.39s 


103.55 g 


1 9+Hs 


2™Bh 


8.80 


8.92 


8.18 


8.93±0.08 


427.56s 


173.58s 


104-76, 


87.98s 


272gi^ 


9.20 


9.61 


8.88 


9.01 


23.83s 


1.49s 


102.36, 


11.8s 




8.60 


8.34 


8.71 


8 80±0 1 


1953 31s 


14779 07s 


102.93, 


1.3min 


264 jjg 


10.20 


10.06 


10.57 




5.94ms 


14.3ms 


10-3.18, 




268Hs 


9.12 


9.59 


9.00 




7.16s 


28 2ms 


10^-24s 




2™Hs 


8.58 


9.44 


8.69 




406.49s 


79.2ms 


lO^-^^s 






9.98 


10.08 


9.20 




22.95ms 


12.1ms 


IOO.6I, 




2^^Mt 


9.00 
10.70 


8.75 
10.13 


9.39 
10.06 


9 30±0 06 

10 02ibl 08 


198.15s 
8 48ms 


1255.21s 
249 99ms 

i.^^ ^ t ^ ^ xxxo 


101.09, 

10-°-4is 


19+°-^2s 
472 6ms 


275i^^ 


10.34 


9.60 


10.06 


9.67 


0.03 s 


3.30s 


10-°-77s 


0.98s 


276j4^ 


10.00 


8.92 


9.93 


9.71 


0.55s 


906 38s 


10-0.05, 


0.72s 


278Mt 


9.30 


8.45 


9.37 


9.55±0.19 


57.68s 


34927.55s 


101.54, 


11.0s 


270Ds 


11.38 


10.41 


10.31 




3.19xl05s 


7.1ms 


10-1.89, 




279Ds 
28lDs 


10.10 
13.34 


8.85 
8.37 


9.68 
8.55 


9.70±0.06 
8.83 


55.67ms 
4.21xl0-^s 


3047.07s 
1.39xl05s 


100.91, 
104.52, 


"•^-0.04^ 

1.6min 


277Rg 


12.04 


10.04 


11.50 




1.31xl0-5s 


0.865s 


10-3-72, 




278Rg 


11.70 


9.89 


11.39 


10.69±0.08 


0.15s 


4.93s 


103.11, 


6.2ms 


279Rg 


11.30 


9.52 


10.92 


10.37 


0.56s 


27.53s 


10-2.36, 


0.17s 
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Continue of Table- in 



Nucleus 


QaCSHF) 


g«(RMF) 


g„(FRDM) 


Q«(exp.) 


T«(SHF) 


r„(RMF) 


TaCFRDM) 


r«(exp.) 


280Rg 


10.80 


9.20 


10.13 


9.75 


19.59s 


557.61s 


100-03s 


5.1s 


28lRg 


14.77 


12.99 


9.37 




1.55xl0-i0s 


1.64x10-^8 


101.90s 




282Rg 


14.60 


8.76 


8.79 


9.00±0.1 


6.29xl0-i0s 


15526.06s 


lO^iOs 


0.748 


282112 


15.87 


9.37 


9.42 




1.04xl0-i2s 


29.46s 


10l-30s 




283112 


14.10 


9.09 


9.01 


9.34±0.06 


6.32xl0-9s 


2666.54s 


1q3.66s 


3.8lo;2s 


284112 


11.48 


10.57 


8.69 


9.17±0.05 


7.15xl0-^s 


0.012s 


103.67g 


9.81^^88 


285112 


8.00 


9.98 


8.59 


9.15 


1.97xl0'^s 


5.44s 


10509s 


298 


294112 


6.12 


7.71 


7.53 




1.69xl0i^s 


2.46x10^8 


108158 




282113 


11.12 


9.89 


9.78 


10.62±0.08 


8.20ms 


21.288 


10l.66s 


88.9m8 


284113 


10.50 


9.68 


9.15 


9.97±0.05 


46.84ms 


94.43s 


IQS.eig 


1.36s 


285113 


5.60 


9.45 


8.97 


9.48±0.11 


5.74xl0i9s 


203.59s 


103-86s 


7.98 


286113 


4.00 


10.34 


8.90 


9.63±0.1 


9.55xl033s 


1.278 


104-44s 


28.38 


286114 


7.60 


10.05 


9.39 


10.19±0.06 


4.57xl08s 


1.268 


10208s 


13+0.01 
"•i-'-o.oi^ 


287114 


8.80 


9.94 


9.31 


9.54 


1.23xl05s 


29.918 


103-39s 


3.88 


288114 


7.70 


8.61 


9.16 


9.84±0.05 


1.71xl08s 


4.83xl04s 


102.80, 


1 Q+3.8 
1-^-0.8^ 


289114 


7.60 


8.66 


8.87 


9.82 


5.33xl09s 


3.67xl05s 


104.86, 


2.6s 


296114 


6.75 


8.32 


8.68 




4.68xl0i2s 


5.864x10^8 


104-43, 




298114 


15.6 


7.91 


8.91 




6.92xl0-i2s 


2.389x10^8 


103.63, 




287115 


7.09 


11.01 


10.25 


10.59 


1.49xl0i2s 


42.8m8 


100.63, 


32ms 


288115 


10.80 


10.75 


10.12 


10.48 


318.88ms 


426.6ms 


101.33, 


250ms 


289115 


10.60 


10.74 


10.03 


10.31±0.09 


482.31ms 


204.84ms 


101.26, 


320ms 


290115 


12.20 


9.23 


9.92 


9.95±0.4 


0.165s 


11050.168 


101.90, 


0.0238 


291115 


8.86 


10.15 


9.66 




87783.45s 


8.458 


102.34, 




290116 


8.90 


11.76 


11.12 


10.84±0.08 


24443.12s 


0.23m8 


10-212s 




291116 


8.50 


10.58 


11.11 


10.74±0.07 


7.77xl06s 


2.25s 


10-1.03, 


18+fms 


292116 


9.40 


9.689 


10.82 


10.71±0.15 


528.65s 


66.091s 


10-1.37, 


33+}fms 


293116 


9.30 


9.679 


10.94 


10.67±0.06 


12923.99s 


825.508 


10-0.59, 


53i^^m8 
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FIG. 5: Same as Fig. 3, but with RMF(NL3*). 
Continue of Table- III 



Nucleus 


Qa(SHF) 


QaCRMF) 


Q^(FRDM) 


Qa(exp.) 


Ta(SHF) 


r„(RMF) 


r„(FRDM) 


T^(exp.) 


293117 


10.70 


9.989 


11.68 


11.03±0.08 


1.07s 


104.71s 




21ms 


294117 


9.40 


10.00 


11.67 


10.81±0.1 


15027.71s 


213.62s 


10-2.05, 


112ms 


297117 


10.74 


9.52 


11.74 




0.844s 


2869.78s 


10-2.55, 




294118 


11.10 


9.60 


12.28 


10.84 


33.36ms 


553.44s 


10-4-24s 


7.1ms 


297118 


4.40 


9.95 


12.10 




5.13xl032s 


555.59s 


10-2.79, 




292120 


12.42 


10.37 


13.89 




9.68x10-^8 


12.30s 


10-6.96, 




293120 


12.24 


11.05 


13.69 




2.81ms 


2.09s 


10-5.51, 




304120 


12.07 


10.19 


13.82 




0.58ms 


42.53s 


10-6.83, 





Here, BE(N,Z) and BE(N-2,Z-2) are the binding energy of the parent and daughter nuclei, respec- 
tively. BE{2, 2) is the binding energy of a particle (^He) i.e. 28.3 MeV OSll . We are using the 
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FIG. 6: Same as Fig. 3, but for 



292,293,304 



120 with RMF(NL3*). 



empirical Viola-Seaborg formula 152 . 



531] for the calculation of Tq which is define as: 



log,o(T„) = + {cZ + d) + /.log, 

V 



(19) 



where Z is the charge number of the parent nucleus and a, h, c, d are the fitting constants defined 
as I52I-I54I1 a = 1.66175, b = -8.5166, c = -0.20228, d = -33.9069. The last term in above equation 
is known as even-odd hindrance factor and is obtained by fits to odd nuclei. It takes the value 

as hlog = 0.0{forZ = even,N = even), hlog = 0.772{forZ = odd,N = even), hlog = 

1.066(/orZ = even, N = odd), hlog = 1.114(/orZ = odd, N = odd). From the above 

1 

equation, one can find that Tq oc 10 v^, hence a small change in creates a large difference in 
Ta- This observation reflects in our results (see Table III). 



IV. SUMMARY 



In summary, we study a broad spectrum of the superheavy nuclei for Z = 105 - 120 in the 
frame work of non-relativistic and relativistic mean field formalisms. The calculations are done in 
an axially deformed coordinates using constant gap BCS pairing. Except Z=\20, all other isotopes 




r(fm) 

FIG. 7: The neutron, proton and total matter density distribution for ^^^117 with RMF(NL3*). 

are synthesized artificially whose properties are studied. From the calculated binding energies, Qa 
and Tq, are estimated. The results are compared with the most acceptable FRDM data as well as 
with the experimental observations, wherever available. Over all comparison of the results are 
quite agreeable with each other in all the three models along with the experimental data. Another 
important observation of our present analysis is the constant value of neutron to proton ratio of 
~ 1.5 for all the observed nuclei. This is enough indication to conclude that the superheavy nuclei 
which will be synthesized in future will have neutron to proton ratio of the order of 1.5. However, 
the binding energy per particle will have much less than the present finding of BE ~ 7 MeV. 

On the basis of the depletion of central density, a bubble structure is studied for some super- 
heavy nuclei in the excited states. This observation is model independent. The maximum depletion 
fraction calculated in the considered nuclei is for ^^^'^^^120. In some of the cases, [Fig.|71 a de- 
pletion of proton density nearer to the center (but not at the center exactly) indicates a special 
type of proton distribution. In this case, the center is slightly bulgy and a considerably depression 
afterward. Again a big hump almost at the mid distance of the center and the surface appears in 
the density distribution for the spherical state. This feature of the nucleus is unique and need more 



investigation. Work in this direction is in progress n55\\ . 



V. ACKNOWLEDGMENTS 

One of the author (MI) acknowledge the hospitality provided by Institute of Physics, 
Bhubaneswar during the work. 



[1] Scharff-Goldhaber, Nucleonics 15, 122 (1957). 

[2] S. G. Nilsson, C. F. Tsang, A. Sobiczewski, Z. Szymanski, S. Wycech, C. Gustafson, I. -L. Lamm, P. 

Moller and B. Nilsson, Nucl. Phys. A 131, 1 (1969). 
[3] G. Herrmann, Nature 280, 543 (1979). 

[4] K. Kumar, Superheavy Elements (TOP Publishing Ltd. 1989); Physica Scripta, T32, 31 (1990); S. 

Hofmann, Rus. Chem. Rev. 78, 1123 (2009). 
[5] W. D. Myers and W. J. Swiatecki, Nucl. Phys. 81, 1 (1966). 
[6] H. W. Meldner, Ark Fys. 36, 593 (1978). 

[7] H. W. Meldner, in Superheavy Elements ed. M. A. K. Lodhi (Oxford: Pergamon) p. 495. 
[8] H. S. Kohler, Nucl. Phys. A 162, 385 (1971). 
[9] S. G. Nilsson, Mat. Fys. Medd. Dan. Vid. Selsk 29, 16 (1955). 
[10] B. R. Mottelson and S. G. Nilsson, Phys. Rep. 99, 1615 (1955); Mat. Fys. Skr. Dan. Vid. Selsk, 8, 1 
(1959). 

[1 1] G. Andersson, S. E. Larsson, G. Leander, S. G. Nilsson, I. Ragnarsson and S. Aberg, Phys. Lett. B 65, 
209 (1976). 

[12] V. M. Strutinsky, Nucl. Phys. A 95, 420 (1967); 122, 1 (1968). 
[13] M. Bhunya and S. K. Patra, 

[14] S. K. Patra, R. K. Gupta and W. Greiner, Mod. Phys. Lett. A 12, 1727 (1997); Tapas Sil, S.K. Patra, 

B.K. Sharma, M. Centelles, and X. Vinas, Phys. Rev. C 69, 044315 (2004). 
[15] Yu. Ts. Oganessian, A. S. Iljinov, A. G. Demin, S. P Tretyakova, Nucl. Phys. A 239, 353 (1975). 
[16] Yu. Ts. Oganessian et al, Phys. Rev. C 74, 044602 (2006). 
[17] Yu. Ts. Oganessian et al, Phys. Rev. C 83, 054315 (2011). 



[18] J. V. Kratz, Heavy Elements and Related New Phenomena, edited W. Greiner and R. K. Gupta (World 

Scientific, Singapore, 1999), Vol. 1, Chap. 4, p. 129. 
[19] A. Sobiczewski and K. Pomorski, Prog, in Part, and Nucl. Phys. 58, 292 (2007). 
[20] J. Decharge, J. -F. Berger, K. Dietrich and M. S. Weiss, Phys. Lett. B 451, 275 (1999); J. Decharge, J. 

-F. Berger, M. Girod and K. Dietrich, Nucl. Phys. A 716, 55 (2003). 
[21] P.-G. Reinhard and H. Flocard, Nucl. Phys. A 584, 467 (1995). 

[22] G. A. Lalazissis, S. Karatzikos, R. Fossion, D. Pena Arteaga, A. V. Afanasjev, and P. Ring, Physics 
Letters B, 671, 36 (2009). 

[23] J.R. Stone, J.C. Miller, R. Koncewicz, PD. Stevenson, and M.R. Strayer, Phys. Rev. C 68 034324 
(2003). 

[24] M. Bender, P -H. Heenen and P -G. Reinhard, Rev. Mod. Phys. 75, 121 (2003). 
[25] M. Dutra, O. Lourenco, J. S. Sa Martins, A. Delfino, J.R. Stone and P. D. Stevenson, Phys. Rev. C 85 
035201 (2012). 

[26] J. R. Stone and P G. Reinhard, Prog. Part. Nucl. Phys. 58, 587 (2007). 

[27] E. Chabanat, P Bonche, P Hansel, J. Mayer, R. Schaeffer, Nucl. Phys.A 627 710 (1997). 

[28] E. Chabanat, P Bonche, P Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A635 231 (1998). 

[29] B.A. Brown, Phys. Rev. C58 220 (1998). 

[30] B. D. Serot, Rep. Prog. Phys. 55, 1855 (1992). 

[31] Y. K. Gambhir, P Ring, and A. Thimet, Ann. Phys. (N.Y.) 198, 132 (1990). 
[32] P Ring, Prog. Part. Nucl. Phys. 37, 193 (1996). 

[33] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

[34] P K. Panda, S. K. Patra, J. Reinhardt, J. A. Maruhn, H. Stocker, and W. Greiner, Int. J. Mod. Phys. E 
6, 307 (1997). 

[35] G. Audi, A.H. Wapstra and C. Thibault, Nucl. Phys. A, 729, 337 (2003). 

[36] P MoUer, J. R. Nix, W. D. Myers, and W. J. Swiatecki, At. Data and Nucl. Data Tables 59, 185 (1995). 
[37] P MoUer, J. R. Nix, and K.-L. Kratz, At. Data and Nucl. Data Tables 66, 131 (1997). 
[38] J. P Maharana, Y. K. Gambhir, J. A. Sheikh and P Ring, Phys. Rev. C 46, Rl 163 (1992). 
[39] S. K. Patra, and C. R. Praharaj, Phys. Rev. C 47, 2978 (1993). 
[40] F. Sarazin et. al., Phys. Rev. Lett. 84, 5062 (2000). 

[41] J. L. Egido, L. M. Robledo, R. R. Rodriguez-Guzman, Phys. Rev. Lett. 93, 282502 (2004). 
[42] J.A.Wheeler (prior to 1950) unpublished. 



[43] H. A. Wilson, Phys. Rev. 69, 538 (1946). 

[44] P. J. Siemens and H. A. Bethe, Phys. Rev Lett. 18, 704(1967). 

[45] C. Y. Wong, Phys. Lett. B 41, 451 (1972). 

[46] K. T. R. Davies, C. Y. Wong and S. J. Krieger, Phys. Lett. B 41, 455 (1972). 

[47] J. Decharge et. al., Nucl. Phys. A 716, (2003). 

[48] M. Grasso et. al., Phys. Rev. C 79, 034318 (2009). 

[49] W. Y. Zao. et. al.. Chin. Phys. Lett, 28, 102101 (2011). 

[50] E. Khan. et. al., Nucl. Phys. A, 800, 37 (2008). 

[51] S. K. Patra and C. R. Praharj, J. Phys. G 23, 939 (1997). 

[52] V.E. Viola, Jr. and G.T. Seaborg, J. Inorg. Nucl. Chem. 28, 741(1966). 

[53] P Moeller, J.R. Nix, and K.L. Kratz, At. Data Nucl. Data Tables 66, 131(1997). 

[54] W.D. Myers and W.J. Swiatecki, Phys. Rev. C 58, 3368(1998). 

[55] S. K. Singh and S. K. Patra (In progress). 



